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The Isolated Photon Cross Section in the Central and Forward Rapidity Regions 

in pp Collisions at y/s = 1.8 TeV 
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A measurement of the cross section for production of single, isolated photons is reported for 
transverse energies in the range of 10-125 GeV, for two regions of pseudorapidity, \n\ < 0.9 and 
1.6 < \rj\ < 2.5. The data represent 12.9 pb _1 of integrated luminosity accumulated in pp collisions 
at y/s = 1.8 TeV and recorded with the D0 detector at the Fermilab Tevatron Collider. The 
background, predominantly from jets which fragment to neutral mesons, is estimated using the 
longitudinal shower shape in the calorimeter. In both pseudorapidity regions the cross section is 
found to agree with the next-to-leading order QCD prediction for 30 < EJ, < 80 GeV. 
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FIG. 1. Distribution of log (-Ei/Utotal) for Ej. = 40 ± 5 GeV central photon candidates (solid points), and the fitted distribu- 
tion (solid curve) made up of Monte Carlo photons (dashes), neutral pions (dots), and rj mesons (dot-dash curve). The neutral 
pion and r\ meson distributions at small log (-Ei/-Etotai) fluctuate due to limited Monte Carlo statistics. 

Direct photon production probes the parton-parton interaction without the ambiguities associated with jet identi- 
fication, fragmentation and energy measurement. In pp collisions at yfs = 1.8 TeV, the dominant mode of production 
for low transverse energy photons is through gluon Compton scattering. The cross section is thus sensitive to the 
gluon distribution in the proton (and antiproton) at low momentum fractions x. With a photon of transverse energy 
> 10 GeV and pseudorapidity r\ < 2.5 (rj = — lntan|, where 9 is the polar angle with respect to the proton beam), 
gluons with momentum fractions as low as x ~ 10~ 3 contribute to the cross section. 

Previous collider experiments, including UA2 |jj and CDF have reported an excess of photons at low E^,(< 
30 GeV) compared with the next-to-leading order (NLO) QCD prediction. This may be explained by additional gluon 
radiation beyond that included in the NLO calculation 0], or by modified parton distributions and fragmentation 
contributions Q. 

This letter presents the first measurement M of the cross section for production of isolated photons in pp collisions 
at y/s — 1.8 TeV with pseudorapidity 1.6 < |?7| < 2.5, and a new measurement (|] of the isolated photon cross section 
in the central region (|?7| < 0.9). 

Photons are identified in the D0 detector j?| using a uranium/liquid argon sampling calorimeter housed in a central 
and two forward cryostats. The calorimeters cover the region of |?7| < 4 and have electromagnetic energy resolution 
(Je/E = 15%/ i/i?(GeV) © 0.3%. In both central and forward regions the electromagnetic section is segmented into 
four longitudinal layers (EM1-EM4) of 2, 2, 7, and 10 radiation lengths respectively; the transverse segmentation is 
into towers of size (in pseudorapidity and azimuthal angle) Ai] x A</> = 0.1 x 0.1 (0.05 x 0.05 at shower maximum 
in EM3). The central and forward drift chambers in front of the calorimeter allow photons to be distinguished from 
electrons and photon conversions by ionization measurement. 
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FIG. 2. Normalized distributions of log (Ei/Etotai) for central electrons from W 
Monte Carlo electrons with Et = 40GeV (histogram). 



ev events (data points), and for single 



The data presented here represent 12.9 ± 0.7 ptT 1 of integrated luminosity recorded during 1992-93. The detector 
used a three-level triggering system. The first level used scintillation counters near the beam pipe to detect an inelastic 
pp interaction. The second level was a hardware trigger which summed the electromagnetic energy in calorimeter 
towers of size A77 x A<j) = 0.2 x 0.2. The data used in this analysis were taken with single tower energy thresholds of 
2.5, 7, and 10 GeV; all except the highest threshold were prescaled. The third level was a software trigger in which 
clusters of calorimeter cells were formed and loose cuts made on shower shape. The cluster energy thresholds used at 
this level were 6, 14, and 30 GeV respectively. 



Photon candidates were selected as follows. Fiducial cuts were applied to select candidate clusters away from the 
calorimeter edges: the clusters were restricted to the regions \r/\ < 0.9 and 1.6 < < 2.5, and in the central region 
were required to be more than 1.6 cm from the azimuthal calorimeter module boundaries. Events where the vertex was 
more than 50 cm from its nominal position were discarded. The resulting acceptance is A — 0.73 ± 0.01 (0.86 ± 0.01) 
in the central (forward) regions. 

The remaining clusters were identified as photon candidates. No drift chamber tracks were allowed in a tracking 
road (A9 x A<fi 0.2 x 0.2 radians) between the calorimeter cluster and the primary vertex. The efficiency of this 
requirement was estimated to be 0.85 ± 0.01 (0.61 ± 0.03) in the central (forward) regions. (The inefficiency is due 
to photon conversions and overlaps with charged tracks from the underlying event.) The photon candidate shower 
was required to have a shape consistent with that of a single electromagnetic shower and to have more than 96% of 
its energy in the electromagnetic section of the calorimeter. The candidates were required to be isolated by a cut 
on the transverse energy in the annular region between 1Z = y 7 Arf + A(f> 2 =0.2 and 7Z = 0.4 around the cluster: 
E^- 0A — E^- ' 2 < 2 GeV. Finally, the missing transverse energy of the event was required to be less than 20 GeV 
to reject electrons from W — > ev decays and events with large amounts of calorimeter noise. The efficiency of these 
last three cuts was estimated as a function of Ej, using a detailed Monte Carlo simulation of the detector and verified 
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FIG. 3. Efficiency-corrected photon purity (V) vs Ej, for central and forward photons. The solid lines show fits of the form 
1 — e ~( a + bE T) anc j t ne dashed lines indicate the range of uncertainty thereon. 



using Z — > ee events. The value obtained was 0.92 ±0.03 (0.77 ±0.06) at £7j = 40GeV for central (forward) photons. 

The primary experimental challenge in the measurement of the direct photon cross section is the extraction of the 
prompt photon signal from the copious backgrounds due to tt° and rj mesons produced in jets which subsequently decay 
to photons. While the bulk of this jet background is rejected by the selection criteria listed above, especially by the 
requirement that the photon candidates be isolated, substantial contamination remains. This comes predominantly 
from fluctuations in the jet fragmentation which lead to a single meson carrying most of the jet energy. If the meson 
has transverse energy above about 10 GeV, the showers from its two decay photons coalesce and mimic a single photon 
shower in the calorimeter. 

The fraction of candidates fulfilling the selection criteria which are genuine direct photons (the purity V) was 
determined using the energy E\ deposited in the first layer (EMI) of the calorimeter. Neutral meson decays produce 
two photons, and so the probability that at least one of them undergoes a conversion to an e + e~ pair in the calorimeter 
cryostat and first absorber plate is roughly twice that for a single photon. Meson showers therefore start earlier than 
photon showers leading to larger E%. A typical distribution of log (Ei/E toia {) is shown in Fig. |l|. The distribution 
is fit as the sum of a photon signal plus 7r° and rj meson backgrounds. Fitting was done separately for the central 
and forward samples for each Ej, bin, using \ 2 minimization, and constraining the fractions of signal and background 
to lie in the range [0,1]. The results presented use a production ratio of ri/ir = 1.0 [||, but all values between 
0.50 and 1.25 give essentially indistinguishable results for V (since the distributions of log (£i/-Etotai) for n° and r) 
mesons are similar). The Monte Carlo calculation combines a detailed simulation of the calorimeter with overlaid 
minimum bias events from data to model noise, pileup and the underlying event. Its ability to correctly model the 
Ei/E to ted distribution has been verified using samples of electrons from W — > ev events taken with the same trigger 
requirements, as shown in Fig. 0. 
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FIG. 4. The inclusive isolated photon cross section a e = d 2 a/dE^,dr/ as a function of photon transverse energy £J, for 
central (circles) and forward regions (squares). The errors are statistical only. The NLO QCD calculated cross sections at, 
using CTEQ2M parton distributions with fj, = Ej,, are shown for comparison. 



The combined statistical and systematic error on the purity V, at each Ej, point (_Ej = 10, 13, 17, 23, 30, 40, 51, 
and 74 GeV), was estimated by inflating by \f)^ the error given by the fit to log(Ea/-Etotai) for that Ej,. The factor 
of \px^ accounts for systematic differences between the Monte Carlo distributions and the data. It was typically 1.3 
in the central region and 1.6 in the forward region. The central and forward photon purities were then corrected by 
the ^-dependent efficiencies and fit to the form: 

V = 1 - e -(«+ 6E ?). (1) 



The data points, fits, and fit errors for V are shown in Fig. g|. 

A second method of purity estimation was used to check the results from the calorimeter energy deposition method. 
It also takes advantage of the difference in conversion probability between single photons and background. In this case 
the material between the interaction point and the central or forward drift chamber (CDC or FDC) is considered as a 
converter and conversions are tagged as tracks with twice minimum ionizing energy using the ionization measurement 
in the CDC or FDC. The results from the two methods are found to be consistent. The conversion method has larger 
statistical errors, since only ~ 10% of photons convert, and therefore it was not used in the fit. 

The differential cross section d 2 cr/dEj, drj is plotted as a function of Ej, in Fig. [|. The next- to-leading order QCD 
calculation was generated using a program due to Baer, Ohnemus, and Owens |9| which includes 7+ jet, 7 + dijet and 
dijet+bremsstrahlung final states. This last is generated from dijet states with collinear bremsstrahlung obtained from 
a phenomcnological photon fragmentation function. In this case, a third collinear jet was created with the remaining 
jet energy, so that the isolation cut could be modeled including jet energy fluctuations. In all cases the parton energies 
were smeared using the measured D0 electromagnetic and jet resolutions. The isolation cut is imposed by computing 



G 



^ 1 

to 

43*0.5 



-0.5 

-1 
1 

0.5 


-0.5 
-1 



. NLO QCD, CTEQ2M, \i = E\ 
\r\ | < 0.9 



, 1o Systematic Error 

IlillllllllliinitiiiNiniiiiiHiuuniNuiiJiiiuiiiuiiiMiiiuillUIIHIIIHIIIIllIlL 



1 .6 < | T| I < 2.5 



□ 



0*. 



-a- 



WIIIIIIIUIIItllL 



1o Systematic Error 



by: 
_ 







_ 



20 40 60 



80 



100 120 

E}(GeV) 



FIG. 5. Difference between the measured isolated photon cross section o e and the NLO QCD prediction a t , normalized to 
the latter. The shaded bands show the magnitude of the combined systematic errors (ler) for each of the two regions. 



the distance 1Z in rj-<p space between the photon and any of the jets, and then rejecting events with a smeared jet 
Et > 2 GeV within 1Z < 0.4 of the photon. (Use of smeared photon and jet energies changes the QCD prediction by 
less than 4% but better represents the actual measurement.) 



The CTEQ2M parton distributions fllCf] and renormalization scale \x 



E T were used. 



If instead the CTEQ2MF 



or CTEQ2MS parton distributions were used, or scales of /j, 
section changes by < 6 %. 



2E T or E T /2 were employed, then the predicted cross 



Figure ^| shows a plot of (a e — <Tt)/&t where tr e and at are respectively the experimental and theoretical values of 
the differential cross section cPa/dEj, drj. The shaded band in the figure shows the magnitude of the systematic error, 
which is estimated by adding in quadrature the uncertainties resulting from the acceptance (~ 1%), the trigger and 
selection efficiencies (3-5%), the photon purity (Fig. 3), the luminosity (5%), and the electromagnetic energy scale of 
the calorimeter (1% in the central, 4% in the forward region). The measured cross sections are in good agreement with 
the NLO QCD prediction for both central and forward regions for moderate transverse energies, 30 < E T < 80 GeV. 
The data points for both the central and forward cross sections lie above the NLO QCD prediction at lower transverse 
energies, but given the magnitude of the systematic error no strong conclusion can be drawn. At the highest transverse 
energies probed, the data for the central region lie below the QCD prediction. Above E T = 74 GeV the photon purity 
V is the result of an extrapolation; the increased systematic error at large Ej, reflects the resulting uncertainty. 

The ratio of forward to central cross sections is shown in Fig. |(| The systematic error on the ratio is estimated by 
adding those for the two regions in quadrature, with the exception of the luminosity uncertainty which cancels. The 
ratio is in good agreement with the NLO QCD expectation for large transverse momenta, but the data lie below NLO 
QCD for xt ;$ 0.04 (E T < 36 GeV). Given the magnitude of the systematic error, no strong conclusion can be drawn. 
This ratio of cross sections, which depends on our unique measurement of the forward photon cross section, could be 
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x T = 2E}/<s 

FIG. 6. Ratio of inclusive isolated photon cross sections d 2 a/dxr dr\ in the forward region to the central region. The NLO 
QCD prediction with CTEQ2M parton distributions is shown for comparison. The shaded band shows the magnitude of the 
combined systematic error. 

used to constrain the gluon distribution function at low x. However, a complete understanding of the origin of the 
low--Ej behavior of the photon cross section is needed before information on the gluon distribution can be extracted. 
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